UNCLASSIFIED 


EOSE 


.  t  o, 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  puxpose 
other  than  In  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  EUid  the  fact  that  the  Govern - 
ment  may  have  formulated,  furnished,  or  In  cmy  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regairded  by  Inqpllcatlon  or  other¬ 
wise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  nay  In  any  way  be  related 
thereto. 


TOPICAL  REPORT 


on 


HEAT  TRANSFER  TO  A  SOLID- 
PROPELLANT  ROCKET-MOTOR  NOZZLE 


to 

AERONAUTICAL  SYSTEMS  DIVISION 
WRIGHT -PATTERSON  AIR  FORCE  BASE 


December  30^  1961 


by 


E.  W.  Ungar 


Contract  No.  AF  33(616)-6358 
Task  No.  73500 


BATTELLE  MEMORIAL  INSTITUTE 
505  King  Avenue 
Columbus  1,  Ohio 


FOREWORD 


This  topical  report  was  prepared  by  Battelle  Memorial  Institute  for 
the  United  States  Air  Force  on  Contract  No.  AF  33(616) -6358.  The  contract 
was  initiated  under  Project  No.  7350  and  Task  No.  73500.  The  work  is 
administered  under  the  direction  of  the  Ceramics  and  Graphite  Branchy 
Metals  and  Ceramics  Laboratory,  Directorate  of  Materials  and  Processes, 
Aeronautical  Systems  Division,  with  Lt.  T.  E.  Lippart  acting  as  Project 
Engineer.  This  report  covers  work  conducted  from  September,  I960, 
through  Jxme,  1961. 
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DEFINITION  OF  SYMBOLS 


A  -  surface  area 

a  -  instantaneous  rate  of  change  of 
b  -  instantaneous  rate  of  change  of  T5 
C  -  constant  of  integration 
Cp  -  specific  heat 

F]^  -  function  defined  by  Equation  (15) 

F2  “  function  defined  by  Equation  (16) 
h  -  heat  transfer  coefficient  defined  by  Equation  (24) 
k  -  thermal  conductivity 

Nu  -Nusselt  number  for  the  gas- side  heat  transfer 
Pr  -  Prandtl  number  for  the  exhaust  products 
q  -  heat  flux 

-  heat  flux  from  the  gas  to  the  exposed  alumina  surface 
(jn  “  heat  flux  from  the  alumina  to  the  nozzle 
r  -  radial  distance 

Re  -  Reynolds  number  of  the  exhaust  products 
T  -  temperature 

Ta  “  temperature  at  the  exposed  alumina  surface 
Tg  -  gas  temperature 
Tg  -  temperature  at  the  nozzle  surface 
t  -  time 

Ug  -  gas  velocity 
Uj  -  liquid  velocity 
V  -  volume  of  ring 

(  Ax)  -  distance  between  centers  of  i  rings 
(Ar)  -  distance  between  centers  of  i  and  j  rings 
a  -  thermal  diffusivity 
6  -  deposit  layer  thickness 
p  -  density 

e  -  T  -  TgO 

Subscripts 

i  -  i^  ring 
j  -  ring 

-  nozzle 

-  initial  time  (t  =  0) 
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HEAT  TRANSFER  TO  A  SO  LID- PROPELLANT 
ROCKET -MOTOR  NOZZLE* 

by 

E.  W.  Ungar 


INTRODUCTION 


The  performance  of  a  nozzle  structure  during  a  solid- propellant  rocket- 
motor  firing  is  intimately  related  to  the  heat  transfer  from  the  hot  exhaust 
products  to  the  nozzle  surface.  Modern  high-energy  solid  propellants  contain 
large  quantities  of  metal  additives.  Thus,  metal  oxides  are  present  in  the 
exhaust,  which  influence  heat-transfer  processes. 

Figure  1  shows  schematically  the  physical  nature  of  the  nozzle  heat- 
transfer  problem.  There  is  an  oxide  layer  of  thickness  6  on  the  nozzle  sur¬ 
face.  In  general,  heat  is  transferred  at  a  rate  qg  to  the  molten  oxide  layer. 

In  addition,  oxide  deposits  from  the  gas  stream  and  also  flows  along  the  sur¬ 
face.  If  the  nozzle  surface  temperature  is  lower  than  the  oxide  melting  point, 
there  is  a  liquid- solid  interface  within  the  oxide  layer.  Heat  is  transferred  at 
a  rate  q^^  from  the  oxide  to  the  nozzle  surface.  In  general,  the  problem  is 
transient,  so  qg  4  qj^,  and  6  is  not  a  constant. 

Early  in  a  rocket-motor  firing  while  the  nozzle  walls  are  cool,  oxide 
accumulates  on  the  walls.  This  oxide  is  cooled  rapidly  by  conduction  to  the 
cold  structure.  After  a  short  while,  the  oxide  layer  grows  to  an  appreciable 
thickness,  thus  providing  some  resistance  to  heat  flow.  When  the  resistance 
to  heat  flow  becomes  sufficiently  large,  it  is  possible  that  the  nozzle  surface 
loses  heat  to  the  bulk  of  the  nozzle  structure  faster  than  it  receives  heat  from 
the  oxide,  resulting  in  a  temporary  drop  in  surface  temperature.  The  hot 
exhaust  products  then  begin  to  reheat  the  oxide.  As  heating  continues,  the 
layer  thickness  begins  to  decrease.  If  the-nozzle  surface  temperature 
reaches  the  oxide's  melting  point,  all  of  the  oxide  on  the  surface  is  molten. 

In  order  to  define  the  exposure  conditions  of  a  nozzle  material,  it 
would  be  desirable  to  define  either  the  nozzle  surface  temperature,  Tg,  and/ 
or  the  heat  flux,  qj^,  as  functions  of  time.  These  quantities  are,  of  course, 
dependent  on  qg,  6,  and  the  physical  properties  of  the  oxide.  To  provide  a 


•  A  supplementary  study  made  in  conjunction  with  Air  Force  Contract  No.  AP  33(616)-6358  which  concerned 
an  investigation  of  nozzle -failure  mechanisms  in  solid -propellant  rockets  for  the  Aeronautical  Systems 
Division,  Wrigfit-Pattenon  Air  Force  Base,  Ohio. 
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FIGURE  1.  SCHEMATIC  OF  PHYSICAL  NATURE  OF  THE 
HEAT-TRANSFER  PROBLEM 
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basis  for  extrapolating  q  to  other  than  the  particular  set  of  test  conditions, 
it  is  desirable  to  relate  u  to  a  heat-transfer  coefficient  and  the  temperature 
difference  between  the  gas  and  the  outer  edge  of  the  o^de  layer. 


SUMMARY 


The  heat  transfer  to  a  small-scale  rocket  nozzle  was  studied  under  the 
conditions  imposed  by  the  combustion  of  two  high-energy  solid  propellants. 
The  exhaust  contain? large  quantities  of  alumi¬ 

num  oxide.  The  heat-transfer  problem  can  be  separated  into  (i^f^he  heat 
transfeiuirom  the  hot  combustion  products  to  a  layer  of  deposited  alumina 
and  ^^the  heat  transfer  from  the  alumina  layer  to  the  underlying  solid 
material,  allpfe  temperature  histories  at  several  locations  in 
nozzles  were  measured  with  high-temperature  thermocouples. 

l-tnupmalUTiaiiiir  thn  — - — ^ 

with  thaWh Baiiia.d -Mhnr  ^ 

j^he  heat-transfer  ra^^rom  the  deposited  alumina  to 
the  nozzles.  Aaaj|KftiK3isaBHMi«t  during  some  periods  of  the  motor  firings^ 
the  heat  transfer  across  the  alumina  layer  could  be  deduced  from  steady- 
state  relations.  The  Uliie^egeig— <  *hs  at  the 


theeat. 


■<^e«wgi4hiB"«etixiimtenl  thiwi 


^e  temper¬ 


ature  at  the  exposed  alumina  surface  and  the  heat-transfer  coefficient  be¬ 
tween  the  combustion  products  and  the  alumina  (MMiitiiee-determined^ijiwjHb^ 

hLUCet'reaga6e3ew»e<<e»wwitgwrtwwiewwwwMMieaeeewi*Mee^4tsi4he»«weiwMeey  'a^-e  ei  with..^ 
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^  turbulent  heat-transfer  theory.  IhiK 


^though  the  alumina-deposit  layer  insulates  the 
nozzle  structure  (and  in  so  doing  strongly  influences  the  tenjperature  history), 
the  effect  of  the  depositing  alumina  droplets  on  gas -side  convective  heat 
transfer  is  negligible.  "Tjy  hniit ■trim afar  deteitminad  in  tln^ 

«tairir;MiMWMappiB«dHta^liii«^«Hip«ta#MWiW#^ozzle  surface-temperature  history  ' 
for  various  ceramic  nozzle  materials.  ^Iculations  successfully  predicted 
the  time  in  the  firing  at  which  melting  or  thermal  degradation  occurred. 


RESEARCH  APPROACH 


Studies  were  carried  out  using  small-scale  solid-propellant  rocket 
motors.  Motor  and  nozzle  configurations  had  been  previously  chosen  for 
materials  evaluations.  Two  high-energy  propellants  were  employed  in  the 
program.  These  were:  (1)  a  composite  propellant  containing  16  per  cent 
aluminxun  (Propellant  A)  and  (Z)  a  double-base  propellant  containing  Z1  per 
cent  aluminum  (Propellant  B). 
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The  over- all  procedure  was  to  measure  temperatures  at  various  posi¬ 
tions  in  a  test  nozzle  during  motor  firings.  The  temperature  histories  were 
used  to  compute  the  heat- transfer  rate  into  the  nozzle  surface,  Chamber- 

pressure  data  were  used  to  determine  the  aluminum  oxide  layer  thickness 
at  the  throat,  6.  Aluminum  oxide  property  data  were  then  used  in  conjunction 
with  a  procedure  described  below  to  compute  the  exposed  alumina  surface  tem¬ 
perature  and  heat  flux,  qg.  A  heat-transfer  coefficient  between  the  gas  and 
the  alumina  layer  was  then  computed  and  compared  with  values  computed  by 
ordinary  methods.  An  inverse  procedure  was  then  used  to  determine  the  tem¬ 
perature  history  of  several  nozzle  materials. 


Calorimeter  Nozzle 


The  short-time  transients  involved  in  a  sol  id- propellant  rocket-motor 
firing  required  that  the  measurement  scheme  be  capable  of  quick  response. 
Surface  temperature  responds  rapidly  to  environmental  changes  and  so  was 
used  to  determine  heat-transfer  rates.  Convergent  test  nozzles  were  fabri¬ 
cated  from  molybdenum  because  its  thermal  properties  are  reasonably  well 
known,  and  because  it  is  machinable,  available,  and  suitable  for  operation 
above  4000  F. 

Figure  2  shows  the  test  nozzle  and  thermocouple  arrangement.  The 
nozzles  were  instrumented  with  ten  tungsten- rhenium  thermocouples.  Each 
thermocouple  wire  had  a  diameter  of  0.  010  inch  and  was  coated  with 
0.  002  inch  of  magnesia  to  provide  insulation  when  placed  in  the  0.  016-inch 
holes  in  the  nozzle.  Junctions  for  five  of  the  thermocouples  were  located 
0.  0125  inch  from  the  exposed  surface  of  the  nozzle,  and  junctions  for  the 
remaining  five  were  0.  0875  inch  from  the  surface.  Each  wire  was  re¬ 
sistance  welded  to  the  bottom  of  a  hole,  two  holes  being  used  for  each 
couple.  The  thermocouple  wires  were  dispersed  around  the  circumference 
®f  the  nozzle  into  16  positions,  each  couple  being  separated  by,  90  degrees. 
Thermocouple  outputs  were  recorded  on  a  high-speed  CEC  galvanometer 
recorder.  The  combined  response  of  the  thermocouple  and  recorder  is  be¬ 
lieved  to  be  better  than  1  millisecond. 

The  placement  of  each  thermocouple  wire  in  a  separate  hole  was  re¬ 
quired  because  of  the  necessity  of  minimizing  temperature- distribution 
disturbances  caused  by  the  presence  of  the  thermocouples.  Because  of  the 
physical  arrangement  of  the  nozzle  on  the  rocket  motor,  it  was  not  possible 
to  drill  the  thermocouple  holes  along  isotherms.  The  disturbance  caused 
by  the  presence  of  thermocouples  was  assumed  to  be  small  because  of  the 
small  size  of  the  holes  and  the  fact  that  both  tungsten  and  rhenium  are  bet¬ 
ter  thermal  conductors  than  molybdenum. 
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figure  2.  MOLYBDENUM  CALORIMETER  NOZZLE  AND 
THERMOCOUPLE  ARRANGEMENT 
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With  the  two-hole  arrangement,  each  thermocouple  consisted  of  a 
tungsten-molybdenum  and  a  molybdenum- rhenixim  junction.  The  molybdenum 
cancels  so  long  as  the  two  junctions  are  close  in  temperature.  The  two 
junctions,  although  at  different  angular  positions,  were  at  the  same  depth 
and  axial  positions.  Measured  temperatures  are,  therefore,  an  average 
of  two  temperatures  measured  at  different  angular  positions.  Thermo¬ 
couples  arranged  in  this  manner  were  calibrated  in  an  induction  furnace  up 
to  ZOOO  F.  The  observed  millivolt  output  was  that  expected  from  an  ordi¬ 
nary  tungsten- rheniiun  thermocouple. 


Numerical  Procedures 


The  determination  of  the  heat  flux  at  the  nozzle  surface,  was 
based  on  a  network  approximation  developed  for  this  study.  In  order  to 
provide  meaningful  results,  the  procedure  had  to  acco’on^’  for  both  the  varia¬ 
tion  in  molybdenvun  properties  with  temperature  and  eucial  conduction  of 
heat  along  the  nozzle. 

The  nozzle  was  considered  as  a  series  of  rings.  A  heat  balance  was 
carried  out  for  each  ring.  The  general  procedure  was  similar  to  those 
described  by  Jakob^^)*  and  Dusinberre'^  with  the  exception  that  the  expres¬ 
sions  for  the  interface  area  between  rings  was  complicated  by  the  nozzle 
shape. 


Figure  3  shows  a  typical  ring.  A  heat  balance  for  the  surface  ring 
can  be  written  as 


q„  iAjAt  +  Ai.i  At 

’  ’  (^^)i-  l,i 

+  Ai  1  kiEi±Jj3!  At  +  Ai  j  At 

’  (Ax)i,i+1  (Ar)i,j 


(1) 


”  Pn  ‘^pn^i  At  "  ^i)  > 

where  the  areas  and  volumes  are  obtained  from  geometric  considerations. 

The  thermocouples  were  placed  at  the  center  of  rings  so  that  the  Tj's  and 
Tj's  could  be  obtained  from  experiment.  Only  the  heat  flux  at  the  surface  was 
sought,  so  the  only  rings  of  interest  were  the  surface  rings  and  the  rings  in  the 
layer  immediately  behind  them.  The  interface  between  the  i  and  j  rings  was 
chosen  so  as  to  be  approximately  parallel  with  the  expected  isotherms.  An 
IBM  650  computer  was  used  for  the  calculation  of  the  local  heat  flux,  from 
the  measured  temperatures. 

*  Numbers  in  parentheses  refer  to  references  listed  at  the  end  of  this  report. 
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The  program  allowed  variable  thermal  properties  to  be  included  in  the 
calculation.  The  thermal  properties  presented  by  Rasor  and  McClelland(^) 
were  used  in  the  computation.  There  is  some  disagreement  between  the 
thermal- conductivity  data  in  Reference  (3)  and  those  presented  in  the  compi¬ 
lation  by  Goldsmith  and  Waterman{4).  Due  to  the  grade  of  the  molybdenum 
used  in  this  program,  the  data  from  Reference  3  were  preferred. 


The  computer  program  described  above  provides  the  local  heat  flux 
at  the  nozzle  surface.  The  aluminum  oxide  layer  on  the  surface  must  be 
analyzed  to  determine  the  gas- side  heat  flux,  qg,  and  the  temperature  at  the 
exposed  alumina  surface,  T^.  The  layer  is  thin  and  can  be  approximated  by 
a  slab  with  an  arbitrary  time  variation  in  temperature  imposed  at  both  faces. 
If,  in  the  temperature  range  of  interest,  the  thermal  conductivity  and  diffu- 
sivity  are  assumed  to  be  constant,  then  the  heat- conduction  equation  becomes 


_  1 

oc  ^ 


(2) 


with  the  boundary  conditions 


r  =  0,  e  =  e^(t)  (3) 

r  =  6,  e  =  eg  (t) 
t  =  0,  e=  e^(r)  , 

where  6  s  T  -  Tg0>  and  r  a  0  is  chosen  to  be  the  gas- side  face  of  the  alumina 
layer. 


The  presence  of  a  possible  solid-liquid  interface  within  the  deposit 
layer  was  neglected  because  of  the  low  latent  heat  of  fusion  of  alumina  and  low 
liquid- film- flow  velocities.  The  layer  thickness,  6,  was  assumed  to  be  con¬ 
stant  for  the  analysis. 


Following  Carslaw  and 

with 

and  at 


Jaeger^^^,  let 


e  =  6j  +  02  > 

(4) 

d^@l  1  d0j 
dr^  a  dt  ' 

(5) 

r  =  0,  01  =  0 

(6) 

o 

II 

CD 

•o 

.  II 
u 

t  =  0,  ©1  =  0o(r). 
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Alto, 


3  02  ^ 

a  S  t 


and  at 


r  =  0,  02  =  ®A(t) 

r  s  6,  @2  =  ®s 

t  =  0,  @2  =  0  . 

Equation  (5)  has  a  solution 

-a  X^t 

01  =  e  (Cl  cos  X  r  +  C2  sin  X r)  ,  ^ 

Introduction  of  the  boundary  conditions  on  r  from  Equation  (6)  leads  to 


6 


1 


wu 

I 


-0(n?Tr^t/6^  •  nTTr 

i  '  am  — • 


n=  1 


At  time  t  =  0, 


00 


?!  a  eQ(r)  =  ^  sin. 


nTir 


6  ’ 


n=  1 


(7) 


(8) 


(9) 


(10) 


(11) 


which  is  simply  a  Fourier  series  for  the  function  0Q(r)  with  constauit  coeffi¬ 
cients  Cjj  given  by 


(r)  sm - dr 

6 


(12) 


Thus, 


01 


^  V  -(xn.^'n^t/5^ 

6  Z/ 

n=  1 


-  .  nTTr  , 

Dq  sin— ^  dr 

0 


(13) 


Thus, 


Duhamel's  Theorem(^)  is  relied  on  for  the  solution  to  Equation  (7). 

62-  e,^<r,  t-T)JdT. 


(14) 


BATTELLE 


MEMORIAL 


INSTITUTE 


10 


Integration  by  parts  yields 


02  =  e^(0)  Fj  (r,  t)  +  08  (0)  F2  (r,  t) 


+ 


t-  T)  +  eg(T)F2  (r,  t-  T)]  dT  , 


(15) 


where  the  dot  denotes  differentiation  with  respect  to  time.  The  function 
F  (r,  t)  represents  the  temperature  at  the  time  t  in  which  the  initial  tempera¬ 
ture  is  zero  and  F^  is  unity  at  r  >  0  and  zero  at  r  a  5,  F2  is  unity  at  r  s  6  and 
zero  at  r  ■  0.  Thus^ 


Fi(r,  t-  T) 


1 


00 


n  s  1 


-on^TT^it  -  T)/6^  .  nTir 
e  ’  sm  — - 

6 


(16) 


t  -  T)  = 


r  2  V  cosnTT  -an^7r2(t- T)/6^ 

6  ^  TT  2/  n  * 
n=  1 


.  nTTr 
sin - 


6 


Combining  Equations  (15)  and  (16)^  leads  to 

00 


n=  1 


2-2^/ A  2 

t/ 6  nTTr 

sin  — ^ 
o 


(17) 


6g(0)  -  0^(0)  +J  (-1)“  Qg-  ©a]  dr} 


Finally, 

*  ■  0  -  f  )  “ A  4  e.  ^  I  f  ,-a.V./42  j-* 


dr  (18) 


n=\ 


■  f  f  ;  X  {  [  "a  (»)  e.  ,0, 


n=  1 


. 4 
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Equation  (18)  is  not  identical  to  the  form  obtained  in  Reference  (5)  but 
is  considered  more  suited  to  this  particular  problem.  The  difference  arises 
because  of  a  difference  in  the  manner  in  which  Duhamel's  Theorem  was 
integrated  (both  forms  are  equivalent).  The  resulting  Equation  (18)  is 
actually  the  limit  of  the  expression  obtained  in  Reference  (5). 

Evaluation  of  the  time  integral  in  Equation  (18)  showed  that  the  prod¬ 
ucts  of  its  value  and  the  value  of  the  negative  exponential  in  time  are  not 
affected  by  0^  or  6g  prior  to  0.  01  second  of  the  point  in  time  under  considera¬ 
tion  for  the  thickness  encountered  and  the  properties  of  alumina.  The  expo¬ 
nential  becomes  quite  small  after  short  times  for  the  thin  alumina  layer  so  that 
the  initial  time  contributions  in  Equation  (18)  are  negligible. 

Equation  (18)  can  be  differentiated  with  respect  to  r  to  yield  the  heat 
flux.  Neglecting  the  initial  time  contributions 


00 

q  =  -k~  =- (0*  -  ej  +  iii  S 

Sr  6  ^  6  4 


cos 


nTTr 


(19) 


n=  1 


X  [e^-  I-  1)”  S,]dT. 


Evaluation  of  Equation  (19)  at  r  =  0  and  r  =  6 ,  and  combination  lead  to 

.t 


qg  =  q^  +  ^  [1  -  (-  1)^1  ^ 


(20) 


n=  1 


For  the  case  of  the  alumina  layer  on  the  nozzle  surface,  evaluation  of 
the  series  shows  that  Equation  (20)  is  closely  approximated  by 


=  9n  +  -  (a  +  b)  pcp6  , 

where  a  and  b  are  the  instantaneous  rates  of  change  of  T»  and  T 

A  I 

The  experimental  data  showed  that 

^n  i  » 


(21) 


so  that  for  the  conditions  in  a  solid- propellant  rocket-motor  nozzle  the 
alumina  layer  could  be  analyzed  on  the  basis  of  a  steady- state  heat- transfer 

equation  with  a  small  correction  of  (a  +  b)  pcp6.  The  temperature  at  the 

exposed  alumina  surface,  T;i^,  could  be  found  from 
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or 
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(22) 


(23) 


These  equations  are  considered  to  be  valid  during  any  time  period  of 
the  motor  firing  that  the  alumina  layer  thickness,  6,  is  approximately 
constant. 


Based  on  a  temperature  difference,  the  heat-transfer  coefficient  is 


h  = 


(24) 


The  calciilation  procedure  was  iterative.  The  value  o/ was  obtained 
from  the  measured  temperatures  by  means  of  the  program  described  earlier. 
The  temperature  was  obtained  from  Equation  (23).  The  alumina 
thermal- conductivity  values  used  in  Equation  (23)  were  those  presented  by 
McQuarrie(6)  and  their  temperature  dependence  had  to  be  included  in  the 
iterative  procedure. 


The  gas  temperature  used  in  the  calculation  should  be  the  adiabatic  wall 
temperature.  This  is  nearly  the  same  as  the  stagnation  temperature,  particu¬ 
larly  in  the  throat  region.  The  temperature  used  in  the  calculation  was  the 
stagnation  temperature  calculated  from  thermodynamic  considerations. 

The  thickness  of  the  alumina  layer  enters  into  Equations  (21)  and  (23). 
Alumina-layer  thicknesses  were  obtained  from  the  chamber-pressure  varia¬ 
tion  in  the  rocket  motor.  Because  of  the  small  size  of  the  nozzle  throat,  the 
chamber  pressure  was  particularly  sensitive  to  the  deposit-layer  thickness  at 
the  throat.  Variation  in  propellant  burning  rate  with  pressure  was  accounted 
for  in  the  calculation.  The  deposit  thickness  computed  for  the  end  of  the  firing 
compared  well  for  each  test  with  the  postfiring  measurement.  Although  de¬ 
posit  thickness  could  be  measured  throughout  the  nozzle,  the  data  were  insuf¬ 
ficient  to  provide  an  alumina-layer  thickness  history  except  at  the  throat. 
Therefore,  the  detailed  heat-transfer-coefficient  analyses  were  carried  out 
only  at  the  throat.  Pressure  dependence  on  alumina  layer  thickness  is  dis¬ 
cussed  further  in  the  Appendix. 


TEST-FIRING  RESULTS 


Five  molybdenum  calorimeter  nozzles  were  fired  for  approximately 
10  seconds  each.  Two  of  these  firings  were  with  Propellant  A  and  three 
were  with  Propellant  B.  A  new  nozzle  was  used  in  each  test  firing. 
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Figures  4  and  5  show  measured  temperatures  typical  of  those  obtained 
with  Propellants  A  and  B,  respectively.  The  thermocouple  numbers  used  in 
Figures  4  and  5  correspond  to  those  in  Figure  2.  The  highest  temperatures 
were  generally  attained  in  the  region  where  the  cone  faired  into  the  cylinder, 
i.  e.  ,  Station  3.  The  temperatures  tended  to  be  only  slightly  lower  at  Sta¬ 
tions  2  and  4.  Axial  conduction  within  the  nozzle  was  important  and  could 
not  be  neglected  in  the  calculation.  Heat  losses  at  the  end  of  the  nozzle 
were  considered  to  be  small  even  when  compared  with  the  eixial  heat  transfer 
within  the  nozzle. 

The  measured  temperatures  increased  rapidly  during  the  tests.  In 
each  of  the  Propellant  B  motor  tests,  the  measured  temperatures  exceeded 
the  4000  F  calibrated  limit  of  the  tungsten- rhenium  thermocouples  at  sev¬ 
eral  locations.  Postfiring  examinations  of  the  nozzles  tested  with  Propel¬ 
lant  B  indicated  that  some  surface  melting  of  molybdenvun  had  occurred. 

The  location  of  the  melting  corresponded  to  the  locations  where  the  thermo¬ 
couples  had  exceeded  their  calibration  limit. 

There  is  a  tendency  for  the  surface  temperature  to  rise  rapidly  early 
in  a  firing,  fall,  and  then  rise  again.  The  second  rise  continues  throughout 
the  remainder  of  the  firing.  This  tendency  was  particularly  pronounced  in 
the  case  of  Propellant  A.  The  temperature  fall  is  probably  associated  with 
the  initial  alumina  deposit,  which  partially  insulAtes  the  surface  layer  while 
it  is  being  cooled  through  heat  absorption  by  the  bulk  of  the  nozzle  which  is 
still  cool. 

Figure  6  shows  a  chamber  pressure-time  curve  for  the  motor  with 
Propellant  A.  The  chamber  pressure  shows  a  minimum  value  following 
ignition,  which  approximately  corresponds  in  time  to  the  temperature  peak. 
Also,  the  chamber  pressure  is  a  maximum  at  approximately  the  same  time 
that  the  temperature  begins  its  final  climb. 

Superimposed  on  the  over- all  trends  are  relatively  high-frequency  tem¬ 
perature  oscillations.  These  are  probably  the  result  of  build-up  and  slough- 
off  of  alumina  from  the  surface.  Similar  oscillations  appear  to  exist  in  the 
measured  pressure,  which  tends  to  confirm  this  explanation. 

The  measured  temperatures  were  used  to  compute  heat- transfer  coef¬ 
ficients  at  the  nozzle  throat  by  the  method  described  previously.  The  high- 
frequency  temperature  oscillations  were  removed  from  the  temperature  curves 
prior  to  the  calculation.  The  calculation  was  carried  out  for  a  time  period 
which  followed  the  alumina-layer  establishment  and  during  which  the  alumina- 
layer  thickness  was  approximately  consteuit.  This  corresponds  to  the  times 
during  which  the  alumina-layer  heat-transfer  analysis  described  previously  is 
valid. 
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FIGURE  4.  TYPICAL  TEMPERATURE  HISTORIES  FOR  PROPELLANT  A 


FIGURE  5.  TYPICAL  TEMPERATURE  HISTORIES  FOR  PROPELLANT  B 


Seconds  After  Ignition 


FIGURE  6.  CHAMBER-PRESSURE  TRACE  FOR 
PROPELLANT  A  MOTOR 
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Table  1  shows  the  mean  value  of  the  computed  heat- transfer  coefficients 
for  Propellants  A  and  B.  These  heat-transfer  coefficients  are  compared  with 
calcvilated  values  based  on  the  method  of  Bartz(^).  Bartz's  approximate  calcu¬ 
lation  procedure  is  based  on  the  similarity  between  nozzle  and  turbulent  pipe 
flows.  In  applying  the  Bartz  method,  no  account  was  taken  of  the  heat  trans¬ 
ferred  to  the  alumina  layer  by  aluminum  oxide  particles.  The  correspondence 
between  the  experimental  and  theoretical  values  of  the  heat- transfer  coefficient 
is  an  indication  that  the  particles  have  no  significant  effect  on  the  heat- 
transfer  process.  Thus  the  heat-transfer  coefficient  should  be  correlated  by 

Nu  =  0.  026  ReO- 8pr0-4.  (25) 

TABLE  1.  HEAT- TRANSFER  COEFFICIENTS 


Propellant 

Computed  From 
Experiment, 

Btu/in.  sec  °F 

Computed  From 
Reference(^), 
Btu/in.  ^  sec  °F 

A 

4.68  X  10-3  ±  8% 

4.  05  X  10-3 

at  700  psig 

at  700  psig 

B 

3.  18  X  10-3  ± 

3.  42  X  10-3 

at  675  psig 

at  675  psig 

Colucci(^)  suggested  a  correlation  similar  to  Equation  (25).  However, 
the  data  presented  by  Colucci  for  an  aluminized  propellant  could  not  be  corre¬ 
lated  by  Equation  (25)  except  in  a  region  of  the  nozzle  exit  cone  where  there 
were  no  deposits. 

The  heat-transfer  coefficient  calculation  was,  of  course,  dependent  on 
the  calculation  of  the  temperature  at  the  exposed  alumina  surface,  T^.  The 
temperature  difference  across  the  alumina  layer  was  found  to  be  of  the  or¬ 
der  of  1000  to  2000  F.  There  was  a  time  period  during  the  firing  when  the 
temperature  T^  was  essentially  constant  and  approximately  equal  to  the 
melting  point  of  alumina.  This  should  be  expected  for  a  portion  of  the  time 
when  the  alumina  thickness  is  decreasing.  Eventually,  the  alumina  surface 
temperature  approaches  the  alumina  melting  point  and  the  temperature 
begins  to  rise  again.  This  occurred  earlier  in  the  firings  with  Propellant  B 
than  with  A  because  of  its  higher  flame  temperature. 
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TEMPERATURE-HISTORY  CALCULATIONS 


The  heat-transfer  coefficients  obtained  from  the  experiments  described 
in  the  previous  section  were  applied  to  the  calculation  of  the  temperature  his¬ 
tory  of  some  ceramic  materials  which  had  been  previously  tested  as  nozzle 
inserts’^  with  the  propellants  used  in  the  present  study.  Since  ceramics 
generally  have  low  thermal  conductivities,  it  is  reasonable  to  apply  a  solution 
to  the  heat- conduction  equation  which  neglects  axial  conduction. 

The  throat  region  of  the  nozzle  is  a  cylinder.  Mayer^^)  showed  that  the 
temperature  in  a  thick-walled  cylinder  heated  on  the  inside  and  insulated  on 
the  outside  is 


where 


00 


n  =  1 


An  = 


-2Na  Ro(Mn) 


(m  n)^R^(/x  n)  -  (/i^+N^)R^(jLt  ) 
'^n  '  n  a'  o'^^n' 


0  =  :^, 
a 


T 

a  2  ' 

a 


U  d 
N  =  — 


R  V-  Y^(Mn) 


(26) 


*  Design  was  the  same  as  that  shown  in  Figure  2. 
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and  a  is  the  cylinder  inside  radius >  b  is  the  outside  radius ^  d  is  the  thickness ^ 
the  J's  are  Bessel  functions  of  the  first  kind^  and  the  Y's  are  Bessel  fvinctions 
of  the  second  kind.  is  an  over-all  heat-transfer  coefficient  which  includes 
both  the  gas -side  heat-transfer  coefficient  h  and  the  thermal  resistaince  of  the 
alumina  layer.  It  is  defined  in  the  usual  manner  by 


U 


o 


1 


where  k^  is  the  thermal  conductivity  of  the  alumina. 

The  eigenvalue  ^  is  determined  by  solving  the  equation 


-Mn 


> 


(27) 


where 


Jl(Mn)  Yi(Mn) 


At  the  inside  surface^  1,  so  that  Equation  (26)  can  be  expressed  as 


00 


.  j-  y  _ 

N,  L  r, 


n  =  1 


O^R|(Mna)-Ro(Mn) 


-  1 


(28) 


Figure  7  is  a  plot  of  Equation  (27)  for  B  =  4  which  corresponds  to  the 
nozzle  used  in  this  study.  It  shows  that;  for  values  of  N^>2.  5,  jLt  is  essen¬ 
tially  independent  of  In  the  rocket  nozzle  vinder  consideration;  Ng^  is  of 

the  order  4. 

<rhe  order  of  magnitude  difference  between  /ij  and  should  also  be 
noted.  A  sample  calculation  showed  that  the  ^J.2,  term  in  the  expansion  of 
Equation  (28)  was  negligible  compared  to  the  /ij  term.  Thus,  only  the  first 
term  in  the  series  of  Equation  (28)  is  required  for  Tg  appreciably  greater 
than  zero.  Additional  terms  are  needed  for  early  times. 

This  analysis  cannot  be  applied  to  the  prediction  of  the  nozzle  surface- 
temperature  history  prior  to  firing  because  the  aliimina-layer  thickness  at 
each  point  in  time  must  be  obtained  from  the  motor  chamber  pressure  by  the 
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FIGURE  7,  EIGENVALUE  VARIATION  WITH  N 
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method  described  in  the  Appendix.  Further,  the  heat-transfer  coefficient 
depends  on  chamber  pressure.  On  the  basis  of  Equation  (25),  the  heat- 
transfer  coefficients  of  Table  1  should  be  corrected  by  an  0.  8  power  for 
variation  in  chamber  pressure.  In  applying  the  results,  however,  a  linear 
chamber-pressure  correction  was  used  to  account  additionally  for  changes 
in  the  burning  conditions. 

The  best  available  thermal-property  data  were  applied  to  the  calcula¬ 
tion.  There  is  considerable  uncertainty  in  the  high-temperature  values  of 
these  properties  for  ceramics. 

Equation  (28)  is  a  solution  to  the  heat-conduction  equation  for  constant 
thermal  properties.  Over  the  temperature  range  of  interest,  thermal  prop¬ 
erties  vary  considerably.  In  order  to  minimize  errors  from  this  source,  the 
values  of  the  thermal  properties  were  averaged  over  the  temperature  range 
encountered  in  the  rocket-motor  firing.  The  mean  values  used  in  the  calcula¬ 
tions  are  shown  in  Table  2. 


TABLE  2.  MEAN  PROPERTY  VALUES  FOR 
VARIOUS  CERAMICS 


Mate  rial 

Density, 
lb/ in.  ^ 

Thermal 

Conductivity, 

10"^  Btu/ sec -in.  -  *F 

Specific 

Heat, 

Btu/ lb  -"F 

Niobium  carbide 

0.  253 

3.  86 

0.  123 

Tantalum  carbide 

0.455 

4.  74 

0.069 

CRB  silicon  carbide 

0.  101 

7.44 

0.  374 

ZT  graphite 

0.074 

7.  80 

0.  430 

The  over-all  heat-transfer  coefficient  was  averaged  at  each  point  in 
time  from  motor  ignition.  Thus,  at  each  point  in  time 

N 

"o,N=|fI  Uo,i  »  (M) 

i  =  1 

where  Uq  ^  is  the  thermal  conductivity  at  each  time  increment  and  N  is  the 
total  number  of  time  increments. 

Time  increments  of  0.  5  second  were  used  in  the  calculation.  An  aver¬ 
age,  taken  in  the  manner  of  Equation  (29),  is  considered  to  be  properly 
weighted  to  account  for  the  variation  in  the  heat-transfer  coefficient  and 
alumina-layer  thickness  during  firing.  The  validity  of  this  approach  improves 
with  time.  It  is  not  considered  valid  during  the  period  of  initial  alumina- 
layer  build-up. 
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Mayer^^)  presented  numerical  results  based  on  a  value  of  S2  =  2  and 
Nj  =  1.  Comparison  with  calculations  of  0^  for  the  nozzles  of  interest  showed 
that  correction  of  the  Mayer  results  by  the  actual  value  of  in  the  manner 
indicated  by  Equation  (28)  yielded  sufficiently  accurate  values  of  6g. 

Figures  8  and  9  show  calculated  surface  temperatures  at  the  nozzle 
throat  for  Propellant  A  and  Propellamt  B,  respectively.  Also  shown  in  the 
figures  are  several  measured  temperature  histories  with  their  fluctuations 
(see  Figures  4  and  5)  removed.  Agreement  between  measured  and  calculated 
temperatures  is  considered  to  be  good  when  thermal-property  differences  (in 
the  case  of  molybdenum)  and  thermocouple  depth  are  accounted  for,  Not  too 
much  faith  should  be  placed  in  the  calculated  temperatures  during  the  first 
2  seconds  of  motor  firing,  i.  e.  ,  during  the  alumina-layer  build-up  period. 
However,  beyond  2  seconds,  the  temperatures  appear  to  be  quite  reasonable. 
The  calculated  temperatures  successfully  predicted  the  time  in  the  firing  at 
which  several  of  the  ceramic  materials  failed  by  melting  or  thermal 
degradation. 


CONCLUSIONS 


The  effects  of  the  aluminum  oxide  deposit  layer  must  be  included  in  the 
analysis  of  nozzle  heat  transfer.  Because  alumina  has  a  relatively  low  ther¬ 
mal  conductivity  and  the  deposit  layer  is  thin,  it  is  possible  to  approximate 
the  heat-transfer  equations  by  the  steady-state  relations  with  a  small  correc¬ 
tion  applied  to  heat  flux. 

Once  the  presence  of  the  alumina  layer  is  properly  accounted  for,  it 
appears  that  a  standard  turbvilent  heat-transfer  correlation  is  reasonably 
correct  at  the  throat.  Detuled  analyses  were  not  carried  out  in  the  con¬ 
verging  section  of  the  nozzle  because  the  deposit-layer  thickness  was  not 
known  as  a  function-of  time.  The  thickness  of  the  layer  depends  on  the  de¬ 
posit  rate  from  the  gas  stream  and  also  on  the  flow  conditions  on  the  nozzle 
surface.  The  alumina-layer  flow  problem  is  considered  particularly  difficult 
because  the  properties  (and  phase)  are  strongly  dependent  on  temperature. 

The  results  of  this  study  were  used  to  successfully  predict  the  throat 
surface  temperature  of  several  ceramic  materieds.  Calculations  based  on 
available  thermal-property  data  successfully  predicted  the  time  in  the  firing 
at  which  melting  or  thermal  degradation  occurred.  For  these  calculations, 
alumina  deposit-layer  thicknesses  had  to  be  inferred  from  measured  pres¬ 
sure  histories.  Completely  a  priori  prediction  of  the  temperature  history 
must  await  an  improved  alumina-layer  theory. 
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FIGURE  8.  CALCULATED  AND  EXPERIMENTAL  SURFACE  TEMPERATURES  OF  NOZZLE  INSERTS 
FIRED  IN  MOTORS  CONTAINING  PROPELLANT  A 
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FIGURE  9.  CALCULATED  AND  EXPERIMENTAL  SURFACE  TEMPERATURES  OF  NOZZLE  INSERTS 
FIRED  IN  MOTORS  CONTAINING  PROPELLANT  B 
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APPENDIX 

PRESSURE  DEPENDENCE  ON  ALUMINA  THICKNESS 


The  mass  flow  of  combustion  products  through  a  rocket  motor  can  be 
approximated  by 

m  =  pAnU  ,  (A-1) 

or  in  terms  of  total  pressure  and  temperature 


m 


Po  A^M 


7+1 

2(7-1) 


(A- 2) 


where  Cj  is  a  constant  which  is  introduced  to  account  for  the  flow  coefficient 
and  for  deviations  from  a  perfect  gas. 


In  addition,  for  solid  propellants,  the  burning  rate  is  proportional  to 
the  chamber  pressure  raised  to  the  exponent  n  (where  for  most  propellants 
n<l).  The  value  of  the  pressure  exponent  is  generally  a  known  property  of 
the  propellant.  Thus,  for  a  constant  burning  area,  which  is  usually  well  ap¬ 
proximated  by  end-burning  propellant  grains, 

m  =  C2  Pq"  .  (A- 3) 

In  order  to  evaluate  the  thickness  of  the  alvimina  layer  on  the  basis  of 
the  measured  chamber  pressure.  Equations  (A-2)  and  (A- 3)  are  combined 
to  eliminate  the  mass  flow  rate.  Thus, 


2 

7rr 


M 


+ 


7+1 

^  2(7-1) 


(A-4) 


It  is  noted  that  at  the  throat  M  =  1.  Also  because  of  the  uncertainties 
in  the  knowledge  of  T^,  7,  and  R,  it  is  convenient  to  consider  terms  contain¬ 
ing  these  quantities  as  part  of  the  constant.  Thus,  Equation  (A-4)  can  be 
rewritten  as 


n-1 


(A- 5) 
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It  is  further  noted  that 


6  =  r^  -  r 


where  r^  is  the  nozzle  throat  radius  prior  to  firing.  Thus, 

n-1 

P  ^ 

6  =  r«  -  -2 -  . 


(A-6) 


(A- 7) 


The  constant  C3  was  determined  on  the  basis  of  the  motoi:- design 
conditions.  This  condition  was  found  to  be  valid  for  firings  with  pyrolytic 
graphite  nozzles.  These  nozzles  had  a  sufficiently  high  surface  temperature 
to  be  free  of  alumina  deposits  during  firing. 

Figure  A-1  shows  plots  of  Equation  (A-7)  for  two  nozzles  which  were 
fired  with  Propellant  A  and  Propellant  respectively.  The  v&lue  of  r 
was  determined  by  measurement  prior  to  firing.  It  should  be  noted  that  the 
sensitivity  of  measured  chamber  pressure  to  alvimina  layer  thickness  is  ap¬ 
proximately  9  X  10 in.  /psi. 
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A-3  and  A-4 


FIGURE  A- 1 .  VARIATION  OF  ALUMINA-LAYER  THICKNESS 
AS  A  FUNCTION  OF  MEASURED  CHAMBER 
PRESSURE 
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